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Abstract 

We study the consequences of introducing individual nonconformity in social inter- 
actions, based on Axelrod's model for the dissemination of culture. A constraint on 
the number of situations in which interaction may take place is introduced in order 
to lift the unavoidable homogeneity present in the final configurations arising in 
Axelrod's related models. The inclusion of this constraint leads to the occurrence of 
complex patterns of intracultural diversity whose statistical properties and spatial 
distribution are characterized by means of the concepts of cultural affinity and cul- 
tural cline. It is found that the relevant quantity that determines the properties of 
intracultural diversity is given by the fraction of cultural features that characterizes 
the cultural nonconformity of individuals. 
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1 Introduction. 



Recently, many dynamical models, inspired by analogies with physical sys- 
tems, have been proposed to describe a variety of phenomena occurring in 
social dynamics [1,2.3,4] . Examples include the emergence of cooperation and 
self-organization, propagation of information and epidemics, opinion forma- 
tion, economic exchanges and evolution of social structures. In this context, 
Axelrod's model [5] for the dissemination of culture among interacting agents 
in a social system has attracted much attention among physicists. 

The concept of culture introduced by Axelrod refers to a set of individual 
features or attributes that are subject to social influence. Agents can interact 
with their neighbors in the system according to the cultural similarities that 
they share. From the point of view of statistical physics, this model is appealing 
because it exhibits a nonequilibrium transition between an ordered final frozen 
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state (a global homogeneous culture) and a disordered (culturally fragmented) 
one [6 , 7.8, 9j . Several extensions of this model have recently been investigated. 
For example, cultural drift has been modeled as noise acting on the frozen 
disordered configurations [TU|. The effects of mass media has been considered 
as external [11] or autonomous [1"2"][1"3"] influences acting on the system. The role 
of the topology of the social network of interactions have also been addressed 
[8lll4|ll5] . Other extensions include the consideration of quantitative instead of 
qualitative values for the cultural traits [16]. These studies have revealed that 
Axelrod's model is robust in the sense that its main properties persists in all 
those cases. In particular, the final frozen states invariably consist of one or 
more homogeneous cultural groups. 

In this paper, we introduce a constraint on the number of situations in which 
interaction may take place, in order to lift the unavoidable homogeneity in 
the final states of the above models. Our model is motivated by the idea 
that generally individuals tend to maintain a minimum degree of identity by 
keeping some cultural features different from those of their neighbors. This 
restriction naturally leads to the persistence of complex patterns of diversity 
in the cultural groups present in the final state of the system. 

The model is explained in Section 2. In Section 3, the results of numerical 
simulations are presented, showing the patterns of diversity in the final frozen 
states. The statistical properties that characterize intracultural diversity are 
calculated in Section 4. Conclusions are presented in Section 5. 



2 Axelrod's model with intracultural diversity. 

Axelrod's model [5] considers a square lattice network of N = L 2 elements 
with open boundaries and nearest neighbor interactions. The state of el- 
ement % is given by a cultural vector of F components (cultural features) 
(s%i, Si2, ■ • • , Sip). Each component Sj/ can adopt any of q integer values (cul- 
tural traits) in the set {l,...,q}. Starting from a random initial state the 
network evolves at each time step following these simple rules: i) An element 
i and one of its four neighbors j is selected at random, ii) If the overlap, 
defined as uj(i,j) = Y.f=i^s if ,s Jf , (number of shared features) is in the range 
< < F, the pair is said to be active with a probability of interac- 

tion equal to u(i,j)/F. iii) In case of interaction, one of the unshared features 
k is selected at random and element i adopts the trait Sjk, thus decreasing in 
one unit the overlap of the pair 

In any finite network the dynamics settles into a frozen state, characterized by 
either u>(i,j) = or = F, Vi,j. Homogeneous or monocultural states 

correspond to = F, and obviously there are q F possible configu- 
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rations of this state. Inhomogeneous or multicultural states consist of two or 
more homogeneous domains interconnected by elements with zero overlap. A 
domain, or a cultural region, is a set of contiguous sites with identical cultural 
traits. Castellano et al. [6] demonstrated that the final states of the system 
experience a transition from ordered states (homogeneous culture) for q < q c 
to disordered states (cultural fragmented) for q > q c , where q c is a critical 
value that depends on F. 

In order to allow for diversity we introduce a parameter Fd such that a pair 
is considered active if the overlap is in the range < uj(i,j) < F — F d , 
with < Fd < F. There is no restriction on which of the Fd features cannot be 
exchanged by an active pair. The case Fd = recovers the original Axelrod's 
model, whereas Fd = F results in frozen configurations for all possible initial 
configurations. The number of possible frozen states is the number of config- 
urations in which neighbors cannot longer interact; thus increasing Fd results 
in an increase of this number. The parameter Fd reduces the number of situa- 
tions in which interactions may take place. In the context of social dynamics, 
the ratio Fd/ F can be seen as a measure of individual nonconformity. 

A cultural region is a set of contiguous sites that possess the same cultural 
vector, whereas a cultural zone is defined as a set of contiguous sites that 
share one or more cultural traits; elements in a cultural zone are said to 
have a "compatible" culture [5]. Cultural zones appear as transient states in 
the original Axelrod's model, but in the final state only cultural regions are 
present. When Fd > 0, cultural zones will usually be present in the final state 
because then contiguous sites have an overlap > F — F d . The model 

can be modified by fixing in advance a subset of features that the elements of 
a cultural zone must share in the final state. 



3 Numerical Results 

As an example of the effects resulting from the inclusion of the parameter 
Fd in Axelrod's model, we shall consider a system of size N = 20 x 20 with 
F = 11, and q = 10, starting from random initial conditions. For Fd = 0, 
the system converges to a homogeneous state, i.e., a single cultural region, 
since q << q c ~ 60. For Fd — 1 the final state consists of a single cultural 
zone possessing few surviving traits. We denote by Kf(s) the number of times 
that the trait value s appears in the fth feature of the cultural vectors in the 
system. That is, Kf(s) = J2iLi^s if ,s- For a particular realization in a system 
of size N = 20 x 20, Table 1 shows Kf(s) for the final state. The first row in 
Table 1 shows that Ki(l) = N and for the remaining traits Kf=\(s ^ 1) = 0; 
that is, all the 400 cultural vectors have reached the value 1 in their first 
feature. Note that all 400 elements also share the traits associated to features 
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/ = 6,9, and 10. In features / = 3,4,7,8, and 11 only 2 different values of 
traits survive; while in features / = 2 and 5 there appear 3 different values of 
traits. The fraction of trait values that disappear during the evolution towards 
the final state of the system is J2f, s ^K f {s),o/(f x (? + 1)). Therefore, for the 
realization in Table 1, about 82% of all the trait values existing in the cultural 
vectors at the beginning have disappeared in the final state of the system. 

In the final state there are 47 different cultural vectors in the system. Unex- 
pectedly, the abundance of these vectors follows a nonuniform distribution. 
Table 2 shows the distribution of the 15 most abundant cultural vectors in 
the final state. These vectors are ranked according to their abundance. The 
abundance of the vector of rank R is denoted by Nr and it is also indicated 
in Table 2. We define the fraction of elements having the cultural vector with 
rank R as Cr = Nr/N. Note that about half of the elements in the final state 
share its cultural vector among the seven most common vectors. 

Figure 1 shows the pattern of the final state of the system. The labels indicate 
the rank of the cultural vector corresponding to each site. Contiguous sites 
having identical cultural vectors are joined by lines. Elements with identical 
cultural vectors tend to form domains, in spite that neighbors with overlap 
F — 1 do not interact. 

The number of surviving traits monotonically decreases during the evolution 
of the system toward its final state. When F d — the number of surviving 
traits in the final state is F if q < q c (all elements are identical and therefore 
there is one trait per feature). As shown above, for F d > the number of 
surviving traits in the final state is greater than F but much smaller than its 
maximum possible value of F x q. On the other hand, the size of the largest 
domain is equal to N for F d = and q < q c , while for F d > the largest group 
involves only a fraction of the elements in the system. Figure 2 shows the 
evolution of both, the number of surviving traits and the fraction of elements 
having the most abundant cultural vector, C\ = N\/N, in a system of size 
iV = 120 x 120, for F d = and F d = 1. 



4 Statistical properties 

Figure 3 shows the average fraction of elements with the most abundant vector 
(Ci) as a function of q, for several values of the parameter F d in a system of 
size N = 40x40 and F = 11. For F d = 0, there exist a critical value q c at which 
the order parameter (C\) exhibits an abrupt transition from a homogeneous, 
monocultural state, characterized by (Ci) = 1, to a disordered, multicultural 
state, where (Ci) C 1 [6|. The value q c is not very sensitive to the variation 
of the parameter F d . For F d > and q < q c , the value of (Ci) is less than 
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1, indicating the presence of intracultural diversity in the single surviving 
cultural zone (there are no neighbors with zero overlap). 

Figure 4 shows the average fraction of cultural vectors (C R ) by rank R, for 
different system sizes. Each curve is the average over 50 realizations of initial 
conditions. For the parameter values used in Fig. 4 (q <C q c ) there is a single 
cultural zone in the final state but several cultural vectors survive having a 
distribution that depends on the system size. The dispersion on each curve is 
larger for low and for high values of R than for intermediate values of the rank. 
The dispersion for low values of R reflects the competition between the more 
abundant vectors during the evolution towards the frozen state. On the other 
hand, the dispersion for large values of R is mostly due to the fluctuations on 
the total number of different vectors present in the final state. 

The frozen patterns in a cultural zone are complex, as shown in Figure 1. 
The cultural diversity can be characterized in terms of the distribution of the 
cultural affinity between any two elements i and j in the system, defined as 
A(i,j) = u(i,j)/F. Figure 5 shows the average distribution of the cultural 
affinities A(i,j) of all pairs of elements in the system, averaged over 10 real- 
izations. The three curves on each panel in Fig. 5 correspond to three different 
values of q. In the case that Fd = 0, the distribution for q < q c consists of 
a single peak at the value A(i,j) = 1, corresponding a homogeneous state, 
while for q > q c a second peak appears at A(i,j) = 0, reflecting the presence 
of multiple domains. As shown in Fig. 5, for F d > intracultural diversity is 
manifested as a wide spectrum of cultural affinities present in the system. As 
Fd increases, the distribution of cultural affinities shifts towards A(i,j) = 0, 
reflecting the increase of intracultural diversity within the cultural zones. 

The spatial distribution of cultural diversity can be characterized by the av- 
erage shell affinity S(r) defined as 

where p(i, r) is the set of elements in a square shell of radius r centered in 
element % (the unit of distance is one site), and n(r) is the number of elements 
on this shell. 

Figure 6 shows S(r) for several values of F d . The average shell affinity S(r) is 
well fitted by the relation 

S(r) ~ S(l) -cdogr. (2) 
We find the scaling S(l) ~ 1 — F d / F and a ~ 1.43-Fd/F for a wide range 
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of values of F, F d and q < q c . The slope 



a 



dS(r) 
dlog(r) 



characterizes the 



gradient of intracultural diversity or cultural dine. Extrapolation of equa- 
tion (2) to S(ro) = allows for a definition of a characteristic distance 



r = exp 



0.7( 



F_ 

Fd 



between two elements that have zero overlap; a con- 
cept that can be applied to the study of the spatial distribution of related 
cultures. These results suggest that the relevant quantity for the description 
of intracultural diversity is the ratio Fd/F. 



5 Conclusions 



In order to allow for individual nonconformity in the Axelrod's model for 
cultural dissemination, we have introduced a parameter Fd that reduces the 
maximum number of shared features for interaction. The inclusion of this pa- 
rameter maintains the main features of the Axelrod's model, corresponding 
to Fd = 0. In particular, there is a nonequilibrium transition from a single 
cultural zone to a multicultural state at about the same critical value q c at 
which the Axelrod's model exhibits a transition from a homogeneous, mono- 
cultural region, to a multicultural state. However, the addition of parameter 
Fd set the stage for the occurrence of complex patterns of intracultural diver- 
sity in cultural zones. The intracultural diversity associated to the constraint 
Fd > is manifested in the appearance of a distribution of the abundance of 
cultural vectors by rank inside cultural zones. We found that this distribution 
is sensitive to the size of the zone, as shown in Figure 4. 

We have introduced the concept of cultural affinity between any two elements 
in order to characterize intracultural diversity in the system. The cultural 
affinity among all the elements in the system for Fd > shows a wide distri- 
bution in contrast to the case Fd = where the cultural affinity can only take 
the values 1 when q < q c , or 1 and when q > q c . 

The cultural affinity between elements separated by a given distance S(r) has 
led us to the concept of cultural cline defined as a = — ■ We found that 
S(r) depends only on the ratio Fd/F since it is well fitted by the relation 
S{r) ~ 1 — F d jF — 1.43F d /Flogr for a wide range of values of F, F d and 
q < q c - 

The model presented here may be useful to describe the emergence of cultural 
gradients such as dialects, gastronomic customs, etc, in geographical areas. 
In the biological context, this model can also be adapted to the study of 
phenotype clines. For instance, it have been proposed that the evolution of 
female mating preferences can greatly amplify large-scale geographic variation 
in male secondary sexual characters and produce widespread reproductive 
isolation with no geographic discontinuity [T?] . 



6 



Acknowledgments 



This work was supported by Consejo de Desarrollo Cientifico, Humamstico 
y Tecnologico of the Universidad de Los Andes, Merida, under grant No. C- 
1275-05-05-B and by FONACIT, Venezuela, under grant No. F-2002000426. 
M.G.C. thanks the Condensed Matter and Statistical Physics Section of the 
Abdus Salam International Centre for Theoretical Physics for support and 
hospitality while part of this work was carried out. 



References 



[I] W. Weidlich, Phys. Rep. 204 (1991) 1. 

[2] S. M. de Oliveira, P. M. C. de Oliveira, and D. Stauffer, Nontraditional 
Applications of Computational Statistical Physics (B. G. Teubner, Stuttgart, 
1999). 

[3] P.W. Anderson, K. Arrow, and D. Pines, The Economy as an Evolving Complex 
System (Addison- Wesley, Redwood, 1998). 

[4] M. San Miguel, V. Eguiluz, R. Toral, and K. Klemm, Computing in Science 
and Engineering 7 (2005) 67. 

[5] R. Axelrod, J. Conflict Res. 41 (1997) 203. 

[6] C. Castellano, M. Marsili, A Vespignani, Phys. Rev. Lett. 85 (2000) 3536. 

[7] D. Vilone, A. Vespignani, and C. Castellano, Eur. Phys. J. B 30 (2002) 299. 

[8] K. Klemm, V. M. Eguiluz, R. Toral, and M. San Miguel, Phys. Rev. E 67 (2003) 
026120. 

[9] K. Klemm, V. Eguiluz, R. Toral, M. San Miguel, Physica A, 327 (2003) 1. 

[10] K. Klemm, V. M. Eguiluz, R. Toral, and M. San Miguel, Phys. Rev. E 67 (2003) 
045101(R). 

[II] J. C. Gonzalez-Avella, M. G. Cosenza, and K. Tucci, Phys. Rev. E 72 (2005) 
065102(R). 

[12] Y. Shibanai, S. Yasuno, and I. Ishiguro, J. Conflict Res. 45 (2001) 80. 

[13] J. C. Gonzalez-Avella, V. M. Eguiluz, M. G. Cosenza, K. Klemm, J.L. Herrera 
and M. San Miguel, Phys. Rev. E 73 (2006) 046119. 

[14] J. Greig, J. of Conflict Resolution 46 (2002) 225. 

[15] M. N. Kuperman, Phys. Rev. E 73 (2006) 046139. 



[16] A. Flache and M. Macy, Los Alamos Arxiv physics/0604201 (2006). 
[17] R. Lande, Evolution 36 (1982) 213. 



7 



TABLE 1 



Kf(s): Traits present in the final state 
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TABLE 2 

More common vectors in the final state 



rank 


occurrence 


accumulated 
fraction 


VECTOR 


1 


58 


0.145 


(1, 


1, 


5, 


2, 


8, 


3, 


o, 


5, 


7, 


5, 


5) 


2 


52 


0.275 


(1, 


1, 


5, 


2, 


8, 


3, 


o, 


1, 


7, 


5, 


5) 


3 


28 


0.345 


(1, 


1, 


5, 


5, 


8, 


3, 


o, 


5, 


7, 


5, 


5) 


4 


23 


0.403 


(1, 


1, 


5, 


5, 


8, 


3, 


1, 


5, 


7, 


5, 


5) 


5 


16 


0.443 


(1, 


2, 


5, 


5, 


8, 


3, 


o, 


5, 


7, 


5, 


5) 


6 


14 


0.478 
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1, 


5, 
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o, 


3, 


1, 


5, 
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5, 


5) 
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13 
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2, 


8, 


3, 


o, 


1, 


7, 


5, 
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9 


12 


0.573 
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1, 


5, 


2, 


4, 


3, 


o, 


1, 


7, 


5, 
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10 


11 


0.600 


(1, 


2, 


5, 


5, 


8, 


3, 


1, 


5, 


7, 


5, 


5) 


11 


11 


0.627 


(1, 


3, 


5, 


2, 


8, 


3, 


o, 


1, 


7, 


5, 


5) 


12 


10 


0.652 


(1, 


1, 


7, 


2, 


8, 


3, 


o, 


5, 


7, 


5, 


5) 


13 


9 


0.675 


(1, 


3, 


5, 


2, 


4, 


3, 


o, 


1, 


7, 


5, 


5) 


14 


8 


0.695 


(1, 


2, 


7, 


5, 


8, 


3, 


1, 


5, 


7, 


5, 


5) 


15 


8 


0.715 


(1, 


2, 


7, 


5, 


8, 


3, 


o, 


5, 


7, 


5, 


5) 
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Fig. 1. Spatial pattern corresponding to a final state for a system of size N = 20 x 20 
with parameters F = 11, = 1, and q = 10. The numbers on each site indicate the 
rank of the cultural vector of that site. The most common vectors (rank R = 1) are 
plotted as solid circles. Sites whose vectors have ranks R = 2, 3, 4 are plotted as open 
circles. In order to facilitate the recognition of patterns a line between contiguous 
sites having identical states are connected by lines. 
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log time 

Fig. 2. Evolution for a system of size N = 120 x 120 with parameters F = 11 and 
q = 10, for Fd = (upper panel) and = 1 (lower panel). The monotonically 
decreasing curve on each panel corresponds to the number of surviving traits (left 
vertical axis) as a function of time. The fluctuating curve on each panel shows the 
size of the largest cultural group N± (right vertical axis) as a function of time. Time 
is measured as the number of iterations that result in a trait exchange. 
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Fig. 3. Average fraction of elements having the most abundant cultural vector, 
(Ci), as a function of q, for several values of F,i in a system of size N = 40 x 40 and 
F = 11. The label on each curve indicates the value of F^. Error bars corresponding 
to one standard deviation are shown on each curve. 
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Fig. 4. Average fraction of cultural vectors (Cr) as a function of the rank R, for 
different system sizes with fixed parameter values = 1, F = 11 and q = 10 < q c . 
Continuous, dashed and dotted lines correspond to N = 60 x 60, 20 x 20 and 
7x7, respectively. Error bars corresponding to one standard deviation are shown 
on several points. 
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Fig. 5. Average distribution of the cultural affinities A(i, j) of all pairs of elements in 
the system, averaged over 10 realizations, for the value of indicated on each panel. 
The continuous, the dashed, and the long-dashed lines on each panel correspond to 
the values of q = 10, q = 62, and q = 66, respectively. System size is N = 80 x 80, 
and F = 11. 
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Fig. 6. The average shell affinity S(r) for different values of Fd as indicated by the 
label besides each curve. System size is N = 80 x 80, F = 11, and q = 10. The slope 
of each curve, a = — ^o^tl) > characterizes the cultural dine. 
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